Algorithms for regional source localization 



Sandra H. Dandach 



Francesco Bullo 



On 

o 
o 



u 

o 



O 

cn 

rn 
O 

ON 

o 



Abstract — In this paper we use the MAP criterion to locate a 
region containing a source. Sensors placed in a field of interest 
divide the latter into smaller regions and take measurements 
that are transmitted over noisy wireless channels. We propose 
implementations of our algorithm that consider complete and 
limited communication among sensors and seek to choose 
the most likely hypothesis. Each hypothesis corresponds to 
the event that a given region contains the source. Corrupted 
measurements are used to calculate conditional posteriors. We 
prove that the algorithms asymptotically find the correct region 
almost surely as long as information is available from three or 
more sensors. We also study the geometric properties of the 
model that make it possible in some situations to detect the 
correct region with a unique sensor. Our simulations confirm 
that the performance of algorithms with complete and limited 
information ameliorates with decreasing noise. 

I. Introduction 

A. Problem description and motivation 

Source localization has assumed increasing interest, and 
has been the subject of study for many researchers. The 
general setting is that one or multiple sources lie in a 
bounded region C, and a group of N sensors divide C 
into N smaller connected localization regions Wi , where 
i € {1,...,N}. They measure a received signal strength 
originating at a source s, the sensors try to cooperatively 
identify the region Wi containing s. We set the problem 
as a multiple hypothesis decision making problem, where 
hypothesis Hi is true if the source lies in the region Wi. 
Maximum a posteriori estimation (MAP) is used as a de- 
cision making tool. We implement the estimation technique 
with an all-to-all and a limited communication algorithm. 
The setting of the problem and the proposed solution prove 
to have some geometric characteristics that we derive later in 
the paper. If properly exploited, these characteristics imply 
the possibility of regional localization with a unique sensor 
for certain source positions. We also prove almost sure (a.s.) 
convergence of both our all-to-all and limited communication 
algorithms. To the best of our knowledge, none of the 
algorithms in the literature provide a similar convergence 
result. 

B. Literature review 

In the classical setting, a number of sensors collaborate to 
locate the exact position of a source. The relation between 
the position of a source and the received signal strength 
(RSS) is described in [1], [2], [3]. RSS indirectly provides 
the distance between the source and a sensor. It is easy to 
formulate a trivial linear algorithm that permits localization 
from the measurements at three sensors. However, such a 
linear algorithm may deliver highly inaccurate estimates of 
the distances, even when the noise is small [4]. On the other 
hand, several authors [5], [6] treat localization as a noncon- 
vex optimization problem. Gradient descent algorithms can 
be used to solve the maximum likelihood estimation prob- 
lems. Other approaches include approximating the nonlinear 



nonconvex optimization problem by a linear and convex one 
and then proposing algorithms for the relaxed problem [7]. 

Following the gradient, and approximating with a linear 
convex problem have limitations. The gradient descent can 
get stuck at local minima far from the correct position, 
leading to the choice of wrong regions, even in the absence 
of noise [8]. Authors in [5] use a method of projection 
onto convex sets. A necessary and sufficient condition for 
the convergence of this algorithm is that the source lies 
inside the convex hull of the sensors. The limitations of these 
methods motivated us to look into regional localization. We 
note that there are instances where the location of the region 
containing a source is all that is needed. A more detailed 
listing of the contributions of this paper is presented below. 

C. Contributions 

This paper presents the source localization problem in a 
setting and formulation that to the best of our knowledge are 
new. We present algorithms based on all-to-all and limited 
communication that require only the computation of integrals 
and therefore present a less computationally exhaustive alter- 
native to the current solutions to the localization problem. We 
also show that as the noise decreases, regional localization 
can be accomplished with a unique sensor for certain source 
positions. We show through an asymptotic analysis that 
choosing Voronoi partitions as localization regions achieves 
zero probability of error in the two sensors case. The most 
important advantage of our formulation is that we are able to 
demonstrate the convergence of our algorithms. We provide 
the proof of a.s. convergence of our algorithms, a step that 
tends to be missing in all of the work presented earlier. 
Finally, the limited communication algorithm is promising 
for the localization problems involving multiple sources. 

D. Paper organization 

The paper proceeds by a problem formulation and an 
explanation of our proposed solution in Section [TT] In Sec- 
tion Un] we derive some asymptotic geometric properties of 
the MAP algorithm when applied to our setting. Section [TVl 
introduces the implementation of the algorithms. The analyti- 
cal proof of almost sure (a.s.) convergence of our algorithms 
is presented in Section [V] Section [Vj] shows our simulation 
results and we conclude in Section IVIII 

II. Problem formulation 

Consider a compact connected environment Ccl 2 . Sup- 
pose there are N disjoint regions Wi, such that U^j^W, = C. 
Suppose also that there are N sensors placed at Xi E Wi 
and that the source located at an unknown location s € 
C, transmits a signal whose power undergoes lognormal 
shadowing described below. 

The average power loss for an arbitrary Transmitter- 
Receiver separation is expressed as a function of distance 
by using a path loss exponent (3 > 2. The power loss 



is proportional to a power of the distance between the 
transmitter and the receiver. For a thorough description of 
signal attenuation models over communication channels, we 
refer the reader to [1], [2], [7]. For reasons to be explained 
shortly, we work with a slight modification of the traditional 
model used in the literature. This model for the received 
power at a sensor i is, P ri — do+ ^-f a _ s \\a , where (3 indicates 
the rate at which the power loss increases with distance, do 
is a nominal distance chosen such that the received power 
in the vicinity of the source is almost equal to P, the 
transmitted power at the source. Note that while this model 
gets rid of the singularity at the source, it converges to the 
same behavior as the classical model used in communication 
literature P ri = p-jSp , when the distance \\xi- s\\ is large. 
Here Pi is the power received at a unit distance from the 
source. Taking noise into account in our model, the received 
power satisfies 



In P.. 



ln(Pdo) - ln(do + \\xi - sf) + m, (1) 



where n, are zero mean, independent and identically dis- 
tributed (i.i.d) white gaussian noise with variance a 2 , each 
associated with a sensor i. The joint probability density 
function of the lnP r = [In P ri , . . . , In P rN ] T , conditioned 
on the hypothesis that the source is at a point y G C is 
given by 



p(lnP ri ,...,lnP r J#) 



1 



(2na 2 ) N / 2 



Ef =1 (InPn-M^ 



exp 



Pd„ 



+ \\xi-y\\ 
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Solving for the exact position of the source requires solving 
for y that will maximize the likelihood of having the received 
observation which becomes the problem of solving for, 



y = arg mm 

v 



N 

E 

i=l 



lnP ri - ln(- 



Pdn 



y\ 



This is a nonlinear nonconvex optimization problem. At- 
tempts to solve it or approximate its solution are a topic 
of great interest. In this paper we look for a regional 
localization, so the conditioning on the exact position y in 
(fJJ is replaced by a regional conditioning. The information 
exchanged between any two communicating sensors are: the 
position of the sensors, the localization regions associated 
with each sensor and the logarithms of the received powers 
(corrupted with log-normal noise). Sensors can share infor- 
mation as soon as they make a measurement. Alternatively 
when the noise level in the communication channel is known 
to be high, it is possible for each sensor to average a set of 
repeated measurements and transmit the averaged logarithm 
of the received power. Averaging helps decrease the noise 
variance, and therefore as we expect and will show later, 
improves the performance. We start by introducing the case 
of one noisy measurement per sensor. 

A. Posterior density with a single noisy measurement 

Since we do not know where the source is, we make a 
worst case assumption on the knowledge of its position s. 



Specifically, we assume that the density of s obeys: 
'1/A, ifseC, 



p(s) 



0, otherwise. 



Here A is the sum of all the areas Aj of Wj, with 
j G {1, . . . , N}. We need to derive the probability density 
conditioned on each hypothesis. 

Lemma II. 1 (Regional conditional density) Let 

z = [lnP ri , . . . , lnP r J T , and note that P(y G Wj) — 
P(Hj) = then 



p(z\y G Wj) 



1 



p(z\y)dy. 



(3) 



Proof: We compute 

d Prob(Z < z,y£ Wj) 



p{z\y G W } 



A 



dz P(y G Wj) 

d S-oo I Wj p( z \y)p(y)dydz 

dz Aj 
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From (0 and (0, we obtain 

p(lnP ri ,...,lnP rN \Hj)P(Hj) = I {2 J 2)N/2 



Ef =1 (InP P| -ln(as 
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dy. (4) 



Similarly, when measurements from only one sensor are 
studied. 

B. Posterior density with aggregated noisy measurements 

In this setting each sensor is allowed to take k repeated 
noisy measurements. Noise independence is assumed be- 
tween sensors and between different samples times for each 
sensor. Then defining 



i=i 



(5) 



the variance of the noise becomes <r 2 (fc) = The regional 
posterior in then becomes 



p{\nV n {k)\Hj)P{Hj) = 



1 



• exp 



A J Wj (2^(fc))i/2 

lnPnW-ln( d0+ |f^||, )) : 
2a 2 (k) 

and the joint conditional regional posterior becomes: 
p(InP ri (fc), . . .,]nr r „(k)\Hi)P(Hi) = ^ [ dy 

A JWi 

11 (w»(fc))v» exp 



dy, 



2a 2 (k) 



Note that, as k — > oo, the noise variance approaches zero, 
and the probability density approaches a delta function. 

Remark II.2 Let 5 be the Dirac delta function. In the 
infinite measurement case, lim^^ a 2 (k) = 0, and the 
probability density satisfies 



p(lnP ri |y) 



lim 



fe^oo (27ra(fc) 2 ) 1 /2 



»(lnP n 



In- 



exp 
Pd 



InP,. 



H- 



Pdg 



2a 2 (k) 



d + \\xi - y\\ 



□ 



Remark II.3 In the sequel, for notational simplicity we will 
treat the aggregated measurement case as if it were identical 
to the single measurement case with the caveat that the 
variance goes to zero. □ 

C. All-to-all information MAP estimation 

In the all-to-all communication (A2A) case, full informa- 
tion is available. Using the conditional probability in (0]i 
MAP selects the hypothesis Hi* according to 



argmaxp(lnP ri , . . . , lnP r „ \Hi)P(Hi 



(6) 



Per Remark III. 31 this selection scheme applies to both the 
single and the aggregated measurement cases. 

Before we proceed to deriving the results in the next 
section, we introduce the definition of the Voronoi diagrams. 

Definition II.4 (Voronoi Diagrams) Given N sensors lo- 
cated at positions {x\, . . . , xn} £ C, we define the Voronoi 
diagram associated with the ith sensor, as follows 

Vi = {x £ C : \\x — Xi\\ < \\x — Xj\\,Vj 7^ i}. 

III. Preliminary properties of regional 

LOCALIZATION FOR ONE AND TWO SENSORS 

In this section we derive certain geometric properties of 
MAP estimation as k — ► oo in (0. These geometric 
properties allow us to conclude the following two results. 
First, for certain source locations, a single sensor suffices to 
asymptotically detect the correct hypothesis. Second, for the 
asymptotic detection problem with two sensors, the selection 
of Voronoi partitions as localization regions leads to exact 
localization. These 2 results should be viewed against the 
fact that, even in the noise-free case, at least 3 non-collinear 
sensors are needed for exact localization. In this section we 
conduct a large sample analysis to prove an interesting geo- 
metric interpretation of the conditional probability densities. 
This analysis recognizes that when k — > oo in (0, the 
Gaussian density approaches a Dirac delta function. Before 
we state the lemma that captures this property, we mention 
a basic property of the Dirac delta function [9]. 

Lemma HLl (On the Dirac delta function) If g : R — * M 

is differentiable and vanishes at positions x 7 , 7 G T, then 



%(*)] - E 



7er 



8{X — Xry) 



(7) 



In keeping with Remark III. 31 consider the situation where 

Pd 



pQnPu\veWj)P(yeWj) 
= A 



5 lnP r 



In 



d + || I* - y\\P 



dy. 



Define the circle C(r,Xi) = {y £ K 2 | \\y - Xi\\ = 
r} and denote its intersection with the region Wj by 
S(Wj,r, Xi) = C(r, Xi) f] Wj. Clearly, this intersection set 
S(Wj,r,X{) is the union of certain arcs of C(r,Xi). Define 
9(Wj,r,Xi) to be the sum of the angles subtended by these 
arcs. Call Hj = (y £ Wj). If we let y = [y\, yj^ and define 

/(yi,y 2 ,Pr i )=lnP ri -lnPd 

+ ln(do + ((x tl - yi f + (x l2 - y 2 ) 2 f /2 ), 

then 



p(lnP n \H j )P(H j ) = ± f S(f( yi ,y 2 ,P ri ))dy 

A JWj 

We are now ready for the following lemma. 



dyi 



Lemma III.2 (The arc-length property) Given a region 
Wj, the conditional probability density satisfies 

pOnP ri \yeW j )P(yeW j ) 

= jJi lnPr >- ln do + Z°-yry y2dm ' 



id, if we let r, = (tt l)' 3 , then 



p(lnP ri |y G Wj)P(y€ Wj 



d 



A(3r. 



,/3-2 



y{Wj,n, Xi ). 



The proof of the lemma is provided in the appendix. This 
lemma can be interpreted as follows. Asymptotically, P ri 
directly provides the circle of radius n centered at a;, 
where the sensor is located. 0(Wj,ri,Xi) is simply the 
angle subtended by the intersection of this circle with Wj. 
The quantity p(\nP ri \y £ Wj)P(y £ Wj) that is used in 
MAP is proportional to this angle. We describe further the 
significance of this result after Lemma HlI.3l In the following 
lemma, we show that having Voronoi partitions as well as 
the MAP estimation algorithm, make the probability of error 
zero in the two sensors case. 

Lemma III.3 (Optimality of Voronoi for 2 sensors) 

Consider two points x\ and x 2 in C C IR 2 . Let V\ and V 2 
be the Voronoi diagrams associated with X\ and x 2 . Take 
s £ C and let r± — \\s — x±\\ and r 2 = \\s — x 2 \\. Then as 
k in © tends to infinity, MAP localization algorithm finds 
the region containing s, with zero probability of error with 
only two sensors. 




Fig. 1. This figure shows two nodes x\ and X2, with a source s £ Vi. 



Proof: The proof of this lemma follows the same 
principle as the proof of Lemma IIII.2I We point out the 
differences below. In figure Q] the source s is indicated 
as one of the two points of intersection of the red and 
green circle. The posterior becomes h(yi , P ri , x x , P r2 , X2) = 
lH(a,w 3 ) 5 (f(a,y2,Pr 1 ,x 1 ))-5(f(a,y 2 ,Pr 2 ,X2))dy 2 , here 
P ri and P,. 2 are the received powers by the sensors. The 
integrals are non zero only at the intersections of two sets 
which we prove to be nothing but the intersection of the 
circles in Figure Q] The definition of Voronoi (Def. III. 4b 
implies in the two sensors case that if s' E C, then s E 
Vi s' G Vl In fact s E V x \\x 1 - s\\ < 

\\X 2 -S\\ <p=> -s'\\ < \\X2 -s'\\ s' E V u if 

s' E C. The other case to consider is s' ^ C. In this case s 
is the unique point of intersection of the circles in C. Both 
ways h is non-zero only in the correct region. ■ 
The two lemmas presented in this section have interesting 
implications. Lemma IIII.2I implies that, for certain source 
locations and as the noise becomes smaller, the MAP esti- 
mation algorithm can determine the correct region containing 
the source with only one sensor. That is true when the 
circle centered at a sensor location with radius j-j is included 
in the region Wj. Lemma IIII.3I on the other hand, gives 
one example where the selection of Voronoi partitions as 
localization regions makes it possible to locate the source 
with only two sensors. This would not have been the case 
with two sensors with general convex regions. 

In general, the noise will not be vanishing and the decision 
needs to be made with a finite number of measurements. 
We will prove in Section [V] that the algorithms presented in 
Section [IV] below converge almost surely as k — > 00 in ©. 

IV. Decision making algorithms with all-to-all 

OR LIMITED INFORMATION 

In this section we present two algorithms based on the 
MAP estimation scheme. The regions can take any shape as 
long as they are compact. In the limited communication case, 
sensors can only talk to their neighbors. We assume there 
is a communication graph that describes the information 
exchange among robots. We consider two cases: the all-to-all 
(A2A) communication case and the limited communication 
case. In the A2A we adopt the complete undirected commu- 
nication graph. In the limited communication case we make 
the following degree assumption: each node has at least two 
neighbors, i.e., each node appears in at least two edges. 

A. All-to-all communication 

In this subsection we present the all-to- 
all communication algorithm (Algorithm #1), 
where we apply the classical MAP estimation 
described in Section III-CI on the complete network. 



Algorithm #1: All-to-all communication MAP 

Network: nodes {1, . . . ,N} with complete communication graph 
State of sensor i is Wi := {In P Ti , Xi, Wi} 
Sensor i executes 

1: For all j € {1, . . . , N}\ {1} 

2: transmit state Wi and receive state Wj 

3: calculate 6j := pQn P ri , . . . , In P rjv \Hj)P(Hj) 

4: find j* := argmaXjg^...^} Oj 

5: Return: decision j* 



Assumption IV. 1 (A2A connectivity and noncollinearity) 

We assume that the graph is complete. That is all nodes can 
communicate with each other. We also assume that at least 
three sensors in the graph are non-collinear. 

B. Limited communication 

In the limited communication algorithm, each sensor ac- 
quires data from its neighbors and calculates a joint condi- 
tional density Each sensor then applies MAP to choose 
the most likely hypothesis. In Algorithm #2 each sensor adds 
the hypothesis of the source being outside its neighborhood, 
computes the corresponding conditional density and com- 
pares it to the densities corresponding to neighboring regions. 
Let TV, = {{j E {1, . . . , N} I j is a neighbor of i} U{i}} 
and let TVj = \Mi | be its cardinality. We described the "source 
outside neighborhood" hypothesis as hypothesis number 
(Hq). We also define lnP rjv , to be the vector composed of 
all measurements In P Ti where i E M- Finally, we mention 
that given an event Hi = x E Wi, the complement event is 
defined by Hi = x E Wf . Computing the density in the 
complement of the neighborhood requires only the addition 
of information about the total region. In fact by applying the 
total probability theorem, we get that 

pQnP rj \Ho)P(H )=pQnP Tj \ (J fli)P( (J Hi) 
= p(lnP ri]yeC )P(y €C)-^2 pQnP rj \HJP(Hi) 

Algorithm #2: Limited communication MAP 

Network: nodes {1, . . . , N} with arbitrary communication graph 
State of sensor i is Wi := {In P Ti , Xi, Wi} 
Sensor i executes 
1: For all j G Mi 

2: transmit state Wi and receive state Wj 
3: calculate := p(lnP rAr . \Hj)P(Hj) 
4: calculate i)O = p(lnP rM *\ H )P(H ) 
5: find j* := argmaXigjv. u{0 } 9i,j 
6: Return: decision j* 



Assumption IV.2 (Limited conn, and non-collinearity) 

We assume that all nodes can communicate to their 
neighbors. We also assume that at least three sensors are 
non-collinear in each neighborhood. 

V. Convergence of algorithms 

In this section we prove that the algorithms presented 
above give the correct decision w.p.l with three or more non- 
collinear sensors. In keeping with Remark [II-3L we examine 
dU as a — > + . We start by stating the following result. 
Lemma V.l (Property of non-collinear sensors) For 
do > and f3 > 0, given a source s E K 2 and three 
non-collinear sensors X\, X2 and x 3 E K 2 , define the 

function f:R^Rby f(z) = £* =1 %+KI? ) 2 
The function f(z) vanishes if and only if z = s. 

Proof: In fact, it is easy to check that the sum is 
zero at z = s. Uniqueness of this solution is verified by 
noting that the sum of the square terms is zero only if all 



the summands are zero. For that to be true one needs to 
find z — (x,y) such that (fx ~ xn) 2 + (y — xa) 2 )^^ 2 = 
- x a ) 2 + (s 2 - x l2 ) 2 f' 2 = (r 2 )?/ 2 , for i G {1, 2, 3}. 
Expanding and subtracting, we obtain that a necessary and 
sufficient condition for uniqueness of the solution is that the 
three points are non-collinear. ■ 
In fact given non-collinear Xi, i G {1, . . . , N} with N > 3, 
a compact region Wj and a source s Wj, there always 
exists 

Dj = min || s - y\\ = dist(Wj,s), 

yeWj 

such that Dj > 0. Also for a compact regions Wj, there 
exists Uj > such that 



Uj > max 

y<=Wj,i£{l,...,N} 



In 



do + \\v-xi\\0 



(8) 



Because of the non-collinearity and the fact that Dj > 0, 
there exists an La > such that 



N 



Lj < min > In 
yeW Z_/ V 

2—1 



|y - ^ 



d 



(9) 



This follows from the lemma above (i.e., from the fact that 
the sum has a unique global minimum at s). Define 



aN 



- U 



j ■ 



(10) 



where a > 1, then we have the following result. 



Lemma V.2 Consider Lj, Uj and r\j as defined in d8), (O 
and ( 1 10b . Suppose the source s is not in region Wj, and 
\ n i\ < Vj f or Q-U i € {1; ■ ■ ■ i -^}- r/ie« 

«-l r . • v (\ d o + 11^-^1^ V nn 
L,- < mm > In— — r + n, . (11) 

a u ^ V do + lls-xj/ 3 7 

Proof: In fact, the sum in (fTTT i satisfies: 

Z— 1 

> /. , - 2UjNr]j - Nrfj = Lj + 2U 2 N 



For simplicity of notation, we choose a = 2 from here 
on. Next, we introduce the final intermediate lemma before 
the main results of this section. 

Lemma V.3 (Upper bound for wrong hypothesis) Given 
Lj, Uj and rjj as defined in ([8}, (O and d 1 Ob , define 



I 3 = p(ln P ri , . . . , In P rN | Hj)P(Hj ) 



1 



/ exp( 



Ef =1 ( 



InR. - In 



A(2ira 2 ) N / 2 

, 2 



Pdn 



do + Wy-XiWf 



2<J 2 



dy. 



If \ni\ < r)j for i G {1, . . . , N}, then for j G {1, . . . , N} 
Aj exp( - Lj/Aa 2 ) 



A(2<ko 2 ) n / 2 
Proof Because of the equality 



lnP r , - In 



Pd Q 



d + \\y-Xi\\P 



In 



d 



d 



the result directly follows from Lemma IV2I and from the fact 
that the surface integral of a function / is upper bounded by 
the surface integral of a constant function g, where g takes 
the maximum value of /. ■ 
We are now ready for the convergence theorem. As usual, 
we define the Q-function Q : K — > R>o by 



Q(x) 



1 



2ir J x 



exp(-y 2 /2)dy. 



Theorem V.4 (Elimination property of wrong hypothesis) 

Consider Xi non-collinear sensors, i G {1, ...,7V} with 
N > 3. Let a be the noise variance. Given a source s (£ Wj, 
then we have 

Prob 

where 

e 3 (a) 



p(ln P ri , . . . , In P rN \Hj)P{H 3 ) < e 3 (a 
Aj exp(-L J / 4cr2 ^ 



> 



Pj{a), 



A(2ira 2 ) N / 2 ' 
Furthermore, as a 



p J {a) = {l-2Q( m /a)) N . 

— > + , we have 

■ 1". 



Proof: From Lemmas IV.2I and IV.3I , we have that 

Prob[p(lnP ri) . . .MP rN \Hj)P{Hj) < ej(a)] 
> Prob [[m, . . -,n N ] T G [-r]j,r]j\ N ] 

N /l 1 
= \\ f 2 _ Prob[n, > r)j] + - - Prob[?i l < -rjj) 

i=i ^ 

= (l-2Q( Vj /a)) N . 

The first inequality comes from the fact that Lemmas IV.2I 
and I V.3 1 hold whenever all |n,| < rjj. The proofs of the two 
limits of ej and \ij are immediate. ■ 
This theorem states that, as a — > + , the probability that the 
joint density function p(lnP ri , . . . ,\nP rN \Hj)P(Hj) takes 
an arbitrarily small value goes arbitrarily close to 1 when Hj 
is not the correct hypothesis. This is so as Q(x) — > as 
x — > oo. To complement the Theorem IV. 41 we prove below 
that for the correct hypothesis, the probability density will 
be lower bounded by a positive term w.p.l. 

Theorem V.5 (Strict positivity for correct hypothesis) 

Consider Xj non-collinear sensors, i G {1, ...,iV} with 
N > 3. Let a be the noise variance. If s G Wj, then we 
have 



Prob[p(lnP r 



MP rN \Hi)P(Hi) > 9(a)] > n(a), 



where 



t/ \ n 7-. i 7-. n A, exp(-L 3 /4cr 2 ) 

*(a)=p(lnP n ,...,lnP r J- £ a(2tt^)^ 



j=l,...,iV 



w= n ^■( <7 )= n (i-zQ^/tr)^. 

3=1, ...,JV 3=1,...,JV 

Furthermore, as a — > + , we /iave 

-v p(lnP ri ,...,lnP rN ) > one/ 0(a) -> 1~. 

Proof: The proof of this theorem follows directly 
from Theorem IV.4I and the total probability theorem. Call 
z = [lnP ri , . . . ,\nP rN ] T . We know from the total proba- 
bility theorem that 

N 



p(z\H i )P(H i )+ Yl K^iW) 



j=l,...,N 

3'¥* 



and, in turn, that 



p(z\H i )P(H i )=p(z)- J2 P(z\H 3 )P(Hj). 



j=l,...,N 



From Theorem I V.4I 



Prob 



piz^Pm^piz)- y, 



j=l,...,N 
3& 



> Prob 



j=l,...,N 



> 



n 



j=l,...,N 



As a — > 0+, *(cr) -> p(z) and 0(ct) -> 1". ■ 
Theorem IV. 5 1 complements Theorem IV.4I in that is shows 
that as a — > + , the probability density conditioned on the 
correct hypothesis is lower bounded by a strictly positive 
term. This event happens asymptotically with probability 1. 

U nder Ass ump tion IIV.1I as MAP follows ©, Theo- 
rems [V3] and [V4] complete the proof of a.s convergence of 
the all-to-all communication MAP algorithm. Similarly for 
the limited communication, under Assumption IIV.21 MAP 
estimation converges almost surely when applied to regional 
localization. 

VI. Simulations 

In this section we show simulation results illustrating the 
type of decision obtained by our algorithms. We also show 
a comparison plot between the all-to-all and the limited 
information algorithms. The plot in Figure |2] shows a correct 




Fig. 2. This figure shows the all-to-all communication algorithm with 25 sensors 
and one source whose region was correctly detected. 



detection of the source. The shaded region corresponds to 
the one detected by the algorithm, the source is shown 
as a star and the sensors as the dots. Figure [3] shows the 




Fig. 3. This figure shows the probability of correct decision making of the all-to-all 
and limited communication cases over 1000 runs, with 20 sensors. The blue curve 
corresponds to the A2A communication algorithm, and the black curve corresponds to 
the limited communication algorithm (with the majority vote). 



results obtained from batches of 1000 runs, with N = 20 
sensors. The plots are obtained as follows: If dj is the index 
of the region chosen by sensor j, we compare it to the 
correct index d*. For the comparison to be fair, and since 
the limited communication algorithm makes decisions only 
about its neighborhood, we compare the two algorithms by 
adding a communication round to the limited communication 
algorithm, where we look at all the sensors that decided that 
the source is in their neighborhood, noting that the decisions 
could be inconsistent, we run a majority vote among the 
aforementioned sensors, and compare the final decision to 
the correct one. The black curve in Figure [3] represents the 
limited communication decisions, while the blue curve repre- 
sents the all-to-all communication decision. Figure [3] shows 
the probability of making a correct decision as a increases. 
It is not surprising that as the noise variance increases the 
probability of correct decision making decreases. 

VII. Conclusion 

In this paper we have presented an all-to-all and a limited 
communication algorithm based on MAP that succeed in 
identifying a region that contains a source. We have also 
presented an asymptotic analysis and derived some geomet- 
ric properties of our algorithms. Those properties had the 
implication that for certain source positions, it is possible to 
solve the regional localization problem with a unique sensor. 
We showed that choosing Voronoi partitions for localization 
regions asymptotically achieves zero probability of error 
in the two sensors case. We were also able to prove that 
when readings from three or more non-collinear sensors are 
available, the algorithms choose the correct region almost 
surely. As an extension to this work we are studying how to 
optimally position the sensors, and choose the localization 
regions, so that we minimize the probability of error of our 
algorithms. 
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Appendix 
Proof: [Proof of Lemma MI. 21 Let 

H(a, Wj) = (j/2 S R such that given y\ = a, [a, y 2 ] € Wj) 



takes the values 



h(yi,P ri ,Xi) 



H(a.Wj) 



5 lnP r . - In 



Pd 



dy 2 



do + \\xi - y\[ 3 , 
5(f(a,y 2 ,P ri ,x i ))dy 2 . (12) 



H{a,W 3 ) 



f(yi,y2,P n ,Xi) = f'(y 1 ,y 2 ,Pr i ,x l ) 



Since 

d 

dy 2 

1 d Q + ((x n -yi) 2 + (x, 2 ~y 2 ) 2 Y 

If we fix yx = a, we can solve for 2/2(0) such that, 
f(a,y 2 (a),P ri ,Xi) = 0. In fact 



fail -yi) 2 + {xi 2 -2/2) 2 ) 2 



f(a,y 2 (a),P ri ,Xi) = 
<^ In P Ti - In ■ 



Pdn 



d a + ((x h - a) 2 + (x i2 - y 2 {a)) 2 ) 2 



^ {x ll - a) 2 + (x i2 - y 2 {a)) 2 = 



pT" 



(13) 



where Ti = \ irp- 1 )c?o J " ■ Observe H(a, Wj) has at most 

two elements satisfying equation (TOT i, one or both of: 



or, 



2/2,1 (a) = x i2 - J rf - {x tl - a) 



2/2,2(0) = x i2 + yjr 2 - (x n - a) 2 , 



(14) 



(15) 



whenever rf > (x^ — a) 2 . Using property (0 of the dirac 
delta function, and substituting with 2/2.1(1) and 2/2,2(0) 
obtained in O and ( TT3T >. (fT2l becomes: 



h(a,P n ,Xi) = I 8(f(a,y 2 ,P n ),Xi)dy 2 , 

! H(a,f) 



( f S{y 2 -y 21 {a)) j 

JH(a,Wj) \f'(a,y 2 ,P r -.Xi)\ U U 2 

if 2/2,1(0) G H(a,Wj) but y 2a (a) $ H(a,Wj) 

f S(y-y 2 ,2(a)) j 

JHia.Wj) \f'{a,y 2 ,P Tz ,Xi)\ U ^ 2 

if 2/2,2(0) G P(a, Wj) but 2/21 (a) £ Pf(a, Wj) 

f iS(2/2-1/2,i(q)) , S(y 2 -y 2t2 (a)) i 

JH(a,W 3 ) \f'{a,y 2 ,P rz ,Xi)\ \f'(a,y 2 ,P rt , Xi )\ a « 2 

if both 2/2,1 (a) and 2/2,2(0) G #(a, Wj) 
Define Ii(a, Wj), the indicator function satisfying 

L(a W ) = l 1 if 2/2,1 (o)ei?(a,W,) 
u ' J ' 1 otherwise 

Similarly define 12(0, Wj), the indicator function satisfying 



Then, (TlZb becomes 

h(a,P ri ,Xi) 



1 



|/'(a,2/2,i(a),P n ,a; i )| 
1 



JiKWj) 



/ 2 (a,Wj). 



|/'(o, 2/2,2(a), Pr,,^) I 

By substituting from ([Pil l, we get 
1 

\f'(a,y 2 {a),P ri ,Xi)\ 
d + ((x h - a) 2 + rf - (x n - a) 2 ) 7 

/V^ ^ (^1 " a ) 2 " «) 2 + *f - (^1 - «) 2 ) 1-1 



d +r; 



/3 



1 d + rf 



/V^"(^-o) 2 rf- 2 /3rf- 2 ^-(^-a) 

Let C,- = {.t e R such that (x, 2/2,1 (a;)) G Wj} and C'j : 
{.t G R such that (x, 2/2,2(2;)) G Wj}. Note that 

igCj^ Jifa:, Wj) = 1 and x G C'j =S> J 2 (a:, Wj) = 1. 

Then, 

V 

lci\ja 



P (hxP ri \y G Wj)P(?/ G Wj) = / h^P^x^dy! 



A J Cj \f'(yi,y2,i(yi),P ri ,Xi)\ 
1 /" 1 



h(yi,P ri ,Xi)dyi+ \ h(yi, P ri ,Xi)dy-i 

i 

dyi 
dyi- 



A Jcj \f'(yi, 2/2,2(2/1), P n ,Xi)\ 



(16) 



Write 



Cj = y A a , with P| A a = , and A a = [a la , a 2 J, 



a = l 



(17) 



= (J A' a , with P| 4, = , and A' a = [a[ a ,a' 2 J. 

a 

(18) 



a=l 



Equation ( [To*] ) can then be written as the sum 

A-p(]nP n \y€W j )P(y€W j ) = 

§ I a. \f'(yi,y2Ayi),Pr„^)\ dyi 
+ § h \r(yi,y2Ayi),Pr t ,x l )\ yi ' 

S flic ^ 

= ]^Ja 2a \f'(yi,y2.Ayi),Pr„^)\ dyi 

+ ^L la \r(yi,y2,2(yi),P^x l )\ dyu 

dp+ff 1 



V do+ ^ • arctan Xil ~ ° — 1° 2 ° 

2^ o j3-2 - arctan 5 rf|o Ia 

Q =i /3< v r r - ( x h -vir 

d Q + rf x,^ - j/i 



+ ° 5 ' • arctan — -. 
Note that 

xa. — a 

arctan ■ 



arctan ■ 



a/ (x^ - a) 2 + (x i2 - 2/ 2 (a)) 2 - - a) 2 

x it - a 7r x l2 - y 2 (a) 

= arctan j— = arctan . 

Xi 2 -y2W 2 x^-a 

The conditional probability density becomes after simplifi- 
cations, 

p(} nP ri \y & W j )P(yeW j ) = ^l.j2 

1 _ arctan **-V2M _1 + arctan ^-w(°iJ 

2 - a 2a 2 - a la 
, d +rf ^ (-k x l2 - y2(a' 2 J n 

\ 7> — 77 y I — eirCtelll ; — 



x l2 - 2/2 K J \ rfo + 



/3 



, a— 1 a—1 



where 9 a and 6*^, are the angles of the arcs in S(Wj, r,, x,) 
described on distinct supports as in ( fTTI ) and ( TT~8T > when 
applicable. ■ 



